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ABSTRACT 


Antenna  coupling  is  of  particular  interest  on  aircraft 
where  the  structure  interacts  with  the  antenna  arrays. 

These  arrays  may  be  extremely  complex  as  new  electronic 
countermeasures  and  signal  intelligence  systems  are  added 
to  aircraft. 

The  thesis  attempts  to  find  a method  to  predict  such 
antenna  interactions.  Four  general  computer  programs  for 
antenna  design  and  radiation  pattern  determination  are 
investigated  ar.d  deficiencies  and  strengths  noted.  These 
programs  primarily  employ  the  moment  method  or  geometrical 
theory  of  diffraction.  Some  interesting  correlations  are 
presented  and  recommendations  on  future  potentially  fruitful 
approaches  are  discussed. 
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I.  INTRODUCTION 


A.  NEED  FOR  STUDY 

Since  the  days  of  Operation  Shoehorn  in  the  early  Viet 
Nam  War  era,  airborne  electronic  warfare  has  experienced  an 
explosive  increase  in  platforms  and  complexity.  All  tactical 
and  strategic  U.S.  aircraft  are  now  equipped  with  various 
electronic  countermeasures  systems  providing  warning  or  pro- 
tection capability.  Likewise,  the  tactical  advantages  of 
near  real-time  signals  intelligence  (sigint)  is  being  recog- 
nized with  projects  such  as  the  S-3  Tactical  Airborne  Signals 
Exploitation  System  (TASES) . 

As  envisioned,  TASES  will  operate  with  the  multi -mission 
carrier  and  provide  passive  early  warning  to  the  local  tactical 
commander.  The  system  will  have  a detection  capability  from 
HF  to  high  J band  and  will  possibly  real  time  data  link  the 
information  directly  to  the  carrier.  Since  the  TASES  aircraft 
will  be  operating  in  a tactical  environment,  it  will  probably 
navigate  and  communicate  using  the  standard  Navy  TACAN,  UHF/ 
ADF,  and  UHF  comm  frequencies.  Also,  there  may  be  additional 
transmissions  from  a data  link  system. 

The  majority  of  these  aircraft  transmissions  will  probably 
be  in  the  UHF  range  (particularly  communication)  and  in  a 
tactical  environment  specific  frequencies  may  change  often 
in  the  course  of  a normal  mission.  To  prevent  saturating 
an  onboard  UHF  detection  system,  some  isolation  must  be  pro- 
vided. The  isolation  can  be  handled  by  selectively  blocking 
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appropriate  frequencies,  judicious  location  of  antennas  for 
minimum  coupling  or  a combination  of  these  techniques. 

The  frequency  blocking  method  requires  the  passage  of 
information  from  the  transmitting  system  to  the  detection 
system  that  a frequency  change  is  occurring.  In  the  simple 
case,  the  pilot  would  tell  the  detection  system  operator  that 
a frequency  change  would  be  made  and  the  operator  would  man- 
ually block  out  the  new  frequency.  This  technique,  while 
inexpensive,  is  fraught  with  potential  human  errors  ranging 
from  forgetfulness  to  misunderstanding.  In  total  degenera- 
tion, the  operator  would  simply  block  out  the  entire  U.S. 
tactical  UHF  range  and  concentrate  on  other  frequencies. 

Going  to  more  complex  frequency  blocking  systems,  the  human 
operators  are  taken  out  of  the  loop,  being  replaced  by  detectors 
and  automatic  switching  systems  with  an  attendant  increase  in 
cost  and  potential  maintenance  problems. 

Clearly,  if  the  systems  could  be  decoupled  at  the  antenna 
environment,  the  frequency  blocking  approach  would  be  unnec- 
essary. However,  in  the  past  there  has  been  no  reliable 
method  of  predicting  antenna  coupling  on  an  aircraft.  Instead, 
antenna  systems  have  been  usually  designed  based  on  past  ex- 
perience and  educated  guesses  and  then  either  tested  on  the 
actual  aircraft  or  on  brass  models.  Often  this  trial  and 
error  method  has  proved  more  costly  than  the  alternative  of 
using  the  complex  blocking  method. 

Ideally,  a computer  simulation  method  can  be  found  that 
will  accurately  predict  the  on  aircraft  performance  and 
coupling  of  antennas. 


10 


B.  BACKGROUND 


1.  Method  of  Moments 

The  performance  of  an  aerial  mounted  on  a conducting 
body  such  as  an  aircraft  depends  on  the  currents  induced  on 
the  whole  structure.  At  frequencies  where  dimensions  of  the 
body  are  comparable  with  wavelength,  the' approximate  current 
distribution  may  be  found  using  the  "method  of  moments"  (MM) . 
This  is  a numerical  technique  for  solving  an  integral  equation, 

in  this  case  the  equation  for  the  current  distribution,  by 

converting  it  to  a matrix  equation  using  a suitable  segmen- 
tation scheme.  Two  segmentation  methods  have  been  used  thus 
far,  namely,  the  patch  and  the  wire  grid  methods  [Ref.  1]. 

In  the  patch  method,  the  body  is  divided  into  a 
number  of  patches,  each  small  in  terms  of  wavelength.  The 
wire  grid  method  uses  a wire  grid  to  represent  the  body  and 
is  based  on  the  principle  that  a wire  mesh  is  equivalent  to 
a conducting  sheet  provided  the  mesh  spacing  is  small  in 
terms  of  wavelength.  The  current  on  an  individual  segment, 
surface  patch  or  wire  segment,  is  assumed  to  be  constant  or 

to  obey  a simple  algebraic  or  trigometric  law.  The  matrix 

equation  is  formed  by  finding  the  reaction  or  mutual  imped- 
ance between  pairs  of  segments,  and  the  excitation  vector 
representing  the  source  points.  The  matrix  equation  is 
solved  to  obtain  the  approximate  current  distribution.  The 
radiation  pattern  may  then  readily  be  obtained  by  summing 
the  contributions  from  individual  segments. 


11 


Applications  of  the  surface  patch  method  to  irregular 
bodies  requires  a surface  representation  with  continuous 
curvature  [Ref.  2)  whereas  straight  wire  segments  are  adequate 
to  define  the  outline  of  the  body,  although  more  segments 
may  be  required.  The  surface  patch  method  also  leads  to 
difficulties  when  applied  to  bodies  with  sharp  discontinuities 
such  as  aircraft  wings.  Thus,  wire  grid  modelling  has  been 
more  extensively  applied  to  aerial  problems  [Refs.  3,  4,  5], 
Knepp  [Ref.  2]  has  used  the  surface  patch  met.iod  to  compute 
the  radiation  pattern  of  an  aerial  mounted  on  a CH-47  heli- 
copter. There  have  been  some  attempts  at  a hybrid  approach 
mixing  the  patch  and  wire  methods. 

Tanner  [Ref.  6]  has  produced  general  computer  methods 
for  solution  of  antenna  problems.  These  methods  include 
direct  integral  approaches,  and  modelling  structures  with 
cylinders  and  wire  mesh  methods.  His  approach  to  the  wire 
mesh  problem  allows  a much  larger  grid  spacing  with  an  at- 
tendant decrease  in  matrix  size  and  computer  time. 

In  all  cases  except  Tanner's,  MM  has  been  considered 
a low  frequency  technique  for  aircraft  limited  to  frequencies 
in  which  the  aircraft  is  of  the  order  of  one  wavelength  in 
size. 

2 . Geometrical  Theory  of  Diffraction 

Another  successful  approach  to  solving  aircraft  type 
problems  is  the  Geometrical  Theory  of  Diffraction  (GTD) . 

GTD  has  been  considered  a high  frequency  technique  and  is 
applicable  to  bodies  that  are  arbitrarily  large  in  an  electrical 
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sense  (Refs.  7,  8,  9] . The  GTD  solution  is  made  up  of  two 
parts,  namely,  wedge  diffraction  and  curved  surface  diffrac- 
tion. 

Using  this  approach,  one  applies  a ray  optics 
technique  to  determine  the  fields  incident  on  the  various 
scattered.  Diffracted  fields  are  found,  using  the  GTD 
solutions,  as  rays  which  are  summed  with  geometrical  optics 
terras  to  give  a far  field  solution.  Rays  from  any  given 
scctterer  will  further  diffract  and  reflect  from  nearby 
structures  causing  various  "higher  order"  terms  to  appear. 
The  yarious  possible  combinations  of  rays  that  interact 
between  the  scatterers  can  be  traced  out  and  determined, 
and  typically  include  only  the  dominant  terms.  Thus  at  high 
frequencies,  only  the  most  basic  structure  features  of  the 
aircraft  need  to  be  modelled. 

3.  Hybrid  Techniques 

There  has  been  effort  recently  to  combine  the  MM  and 
GTD  methods.  This  has  included  approaching  the  problem  by 
expanding  the  MM  technique  to  include  GTD  methods  [Ref.  Sj 
and  the  reverse  approach  of  extending  GTD  via  MM  by  Burnside 
at  Ohio  State  University. 

Both  of  these  hybrid  techniques  allow  investigation 
of  problems  that  were  not  solvable  with  MM  or  GTD  alone. 

The  results  have  just  been  published  recently  and  appear 
encouraging . 

A further  discussion  of  all  of  these  techniques  will 
follow  in  the  analysis  section.  Each  has  its  capabilities 
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and  limitations  and  will  be  examined  in  light  of  the  UHF 
coupling  problem. 

C.  STATEMENT  OF  THE  PROBLEM 

The  remainder  of  this  paper  examines  the  various  computer 
programs  available  to  model  the  coupling  of  two  UHF  stub 
antennas  on  an  aircraft.  The  UHF  antennas  are  modelled  as 
nine  inch  (.23  meter)  whips  with  a bottom  feed.  It  is  recog- 
nized that  actual  UHF  blade  antennas  are  multi-element  and 
driven  differently  (for  broad  band)  but  consultation  with 
NADC  (telephone  conversation  with  Joe  Miller)  indicates  that 
it  is  common  to  model  the  UHF  antenna  this  way.  Also,  NADC 
indicated  that  a realistic  termination  impedance  would  be 
seventy  ohms  and  the  gain  of  standard  UHF  antennas  is  about 
1.8  db  ove  isotropic. 

Three  hundred  megahe  tz  is  the  only  frequency  examined. 
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II.  METHOD  OF  ANALYSIS 


A.  DEFINITION  AND  THEORY  OF  COUPLING 


In  regard  to  antenna  problems,  coupling  is  defined  as  an 
indication  of  the  quality  and  quantity  of  energy  received  at 
one  antenna  from  another  antenna.  It  is  conveniently  expresse-’. 
in  a coupling  coefficient  defined  as 


C = 


*7 


(1) 


where  Pr  is  the  total  received  power  and  Pt  is  the  total 
transmitted  power.  The  coupling  coefficient  relates,  for  a 
given  radiator  aid  receptor,  the  power  received  to  the  power 
transmitted.  The  coupling  coefficient  is  a function  of  the 
source,  path  and  receptor  parameters. 

For  an  ideal  case,  conventional  calculation  techniques 
can  be  used  to  determine  the  far  field  coupling  coefficient. 
Looking  at  coupling  between  a given  set  of  antennas,  the 
power  density  at  any  point  in  space  due  to  the  transmitting 
antenna  can  be  expressed  as 

p B Pt  Gt(e,4>,*) 

d 4 it  Rz  (2) 


where  Gt  describes  the  transmitting  antenna  gain  as  a func- 
tion of  frequency  and  spatial  position.  The  effective  capture 
area  of  the  receiver  antenna  is 

Gr(0,4>,A)A2 

Ar  " 

4 it  (3) 
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where  Gr  describes  the  gain  characteristics  of  the  antenna. 
The  total  received  power  is 


P ■ PtGt(8,$,A)Gr(v' 
r (4  ir  R)  2 


(4) 


And  the  coupling  coefficient  is 

c . Fr  = GtC9,^,l(Gr(6,4>>A)X2 

*1  (4  it  R)  2 (5) 

Note  that  the  coupling  coefficient  is  a function  of  the 

product  of  the  two  antenna  gains. 

Consider  two  short  dipoles  (length  less  than  a quarter 

wavelength),  parallel  in  space.  For  this  case 

Gt(e.$,A)  • Gr(0,$,X)  = 3/2  sin2  q (6) 

so  that  the  coupling  coefficient  becomes 

C . (3/2)  2 (sin  er>2 

(4  a R)  2 (7) 

It  is  interesting  to  note  that  the  NADC  rule  of  thumb 
for  blade  antennas  of  1.8  db  gain  is  virtually  identical 
to  the  3/2  coefficient  in  equation  (6) ! 


C « vm  T)  «■*  <-»  1/M  rnflf  1 hi  4-  U A A 4 A A «•  A A A £ MANAMA  1 A 

i / wih  iwvaWo;  [ i\c  i . iuj  , tac  cucvuvc  aiua  vi  a uiuuupui  c 


on  a ground  plane  is 
>r 


A. 


Gr  A2 


16  » 


(8) 


hus  it  can  be  seen  that  equation  (7)  also  applies  to  a 
short  monopole  in  a ground  plane. 

For  a one-quarter  wavelength  monopole  on  a ground  plane 
the  gain  is 


3.28 


cos(a/2  cose) 
sin  6 


-*  2 


(9) 
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and  the  coupling  coefficient  becomes 


C 


(3.28) 


[£ 


cos(tt/2  cos  9) 
sin  9 


* 

A2 


2 2 ( 4 it  R) 


(10) 


B.  METHOD  OF  MOMENTS  (AMP)  GROUND  PLANE  ANALYSIS 

MB  Associates  Information  Systems  Division  of  California 
has  developed  a comprehensive  program  using  wire  grid  modelling 
known  as  AMP  (Antenna  Modelling  Program)  (Ref.  11]  which  was 
available  for  use  at  the  Naval  Postgraduate  School.  This 
program  was  developed  under  a government  contract  as  a general 
aerial  design  tool. 

The  AMP  program  is  based  on  a rigorous  integral  equation 
for  the  antenna  current  which  is  solved  numerically  using 
matrix  manipulation  techniques.  A sinusoidal  interpolation 
scheme  is  used  to  define  the  current  distribution  on  indi- 
vidual wire  segments.  The  current  is  made  up  of  three  terms. 

Ij  the  current  on  the  jth  segment  is  given  by 


It  (s)  = A-j  + Bi  sin  k (s  - Si)  + Cj  cos  k (s  - Si) 

) ) (11) 

where  s is  the  distance  measured  along  the  wire,  for  s in 
sj  defines  the  segment  and  sj  is  the  midpoint  of  the 
jth  segment,  k = 2tt/A  , A is  the  free  space  wavelength.  The 
coefficients  (A,  B,  and  C)  are  obtained  by  enforcing  the 
requirement  that  the  extrapolated  current  from  a given  seg- 
ment matches  the  midpoint  current  in  two  adjacent  segments 
(Ref.  3] . The  average  midpoint  current  is  used  for  multi- 
segment junctions.  The  integration  for  the  sine  and  cosine 
terms  in  the  impedance  calculation  have  analytic  solutions 
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but  numeric  integration  is  required  tor  the  constant  current 
term.  Kirchhoff's  current  law  is  enforced  at  all  intersec- 
tions. 

The  sinusoidal  interpolation  technique  gives  good  results 
for  wire  aerials  where  segments  of  equal  length  can  be  used 
and  where  there  are  no  multi-segment  junctions  because  the 
current  discontinuities  are  then  very  small.  However,  where 
there  are  junctions  between  several  segments  of  different 
length  there  is  a larger  discontinuity.  In  modelling  irreg- 
ular bodies  this  must  be  kept  in  mind. 

This  approach  is  valid  for  any  frequency  and  any  structure 
size  to  wavelength  ratio  [Ref.  11) . However,  it  is  recommended 
to  keep  segment  lengths  to  less  than  one-tenth  wavelength  and 
grid  spacings  for  solid  surfaces  to  about  the  same  size. 
Reference  11  leavts  open  the  prospect,  however,  that  somewhat 
longer  segments  may  be  acceptable  on  long  wires  with  no  abrupt 
changes.  It  should  be  remembered  that  the  size  of  the  segments 
determines  the  resolution  in  solving  for  the  current  on  the 
model  since  the  current  s computed  at  the  center  point  of 
each  segment. 

With  the  above  restrictions  in  mind,  the  AMP  program  has 
most  commonly  been  applied  to  wire  antennas  and  problems  in 
the  HF  and  UHF  bands. 

AMP  offers  many  options  in  modelling  the  antenna  environ- 
ment including  placing  other  wire  structures  nearby,  using 
series  or  parallel  R-L-C  circuits,  terminal  networks  and 
transmission  lines. 
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As  a check  on  the  theoretical  results  of  coupling 
coefficients  in  the  previous  section,  a series  of  AMP  runs 
were  made  usir*  .23  meter  monopoles  of  .00C1  meter  diameter 
on  a ground  plane.  The  antennas  were  split  into  15  equal 
segments.  The  transmitting  antenna  was  fed  at  the  base  with 
a voltage  source  and  the  receiving  antenna  was  terminated  at 
the  base  with  a 70  ohm  resistance.  The  antennas  were  spaced 
at  various  distances  and  the  transmitted  and  received  powers 
were  calculated. 

The  AMP  program  gives  a direct  output  of  radiated  power. 
Received  power  was  calculated  from  the  output  of  the  current 
in  the  70  ohm  element  from  Pr  * i?  R (R  is  resistance,  i.e., 
70  ohm,  here).  The  coupling  coefficients  were  calculated 
with  equation  (1) . 

See  Table  I for  a comparison  of  these  results.  The 
comparison  for  ground  planes  show  reasonable  agreement  with 
theoretical  results.  The  lower  coupling  of  the  AMP  antennas 
is  probably  due  to  the  less  than  optimum  impedance  match  on 
the  receiver  antenna.  Th5  constant  ratio  of  error  in  the 
AMP- divided -by- the -quart er-wave- theoret ical-result  column 
lends  support  to  the  hypothesis  that  the  difference  is  con- 
sistant  and  due  to  the  difference  in  the  model.  At  the  one 
meter  point,  the  near  field  solution  included  in  AMP  can  also 
be  seen  in  the  error  ratio  change. 

The  results  for  the  ground  plane  analysis  should  be 
applicable  to  antenna  arrays  on  large  flat  surfaces  where 
the  antennas  are  mounted  at  least  2 to  3 wavelengths  from 
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an  edge.  The  upper  surface  of  a C-5A  wing  or  its  lower 
fuselage  arc  examples  of  such  cases. 


Distance 

Meters 


1.0 

2.7 

5.0 

8.0 
10.0 
12.2 


TABLE 

Theoretical  Coupling 
Compared  to  AMP 


Ceqn  (7)* 
Short 
Monopole 
Ground 
Plane 

Ceqn  (10) 
/ 4 

Monopole 

Ground 

Plane 

-18.46 

-17.67 

-27.09 

-26.29 

-32.44 

-31.65 

-36. S2 

-35.73 

-38.46 

-37.67 

-40.19 

-39.39 

Coefficients 

Results 


AMP 
. 23m 

Monopole 

Ground 

Plane 

AMP/C  (10) 

-19.32 

-1.65 

-27.51 

-1.22 

-32.91 

-1.26 

-36.97 

-1.24 

-38.91 

-1.24 

-40.63 

-1.24 

•Coupling  coefficients  expressed  in  db. 


An  attempt  was  made  to  use  AMP  to  calculate  coupling  over 
a plate.  The  220  segment,  one  meter  plate  described  in  section 
II.  E.  was  used  with  a source  antenna  at  its  center  and  a 
receive  antenna  at  various  locations  on  the  plate.  The  source 
and  receive  antennas  were  the  same  .23  meter  whips  used  in  the 
ground  plane  analysis  above.  Comparative  ground  plane  and 
theoretical  results  using  equation  (10)  were  calculated.  The 
results  are  shown  in  Table  II. 

TABLE  II 

Coupling  Coefficients  over  a Grid 


Distance 

Meters 

C • 

Eqn  (10) 

C 

AMP  g.p. 

AMP  grid  Cgrid/CR.P 

.25 

-S.63 

-9.54 

-43.99  (on  plate) 

-34.35 

.50 

-11.65 

-14.26 

-50.47  (edge) 

-36.21 

.707 

-14.66 

-16.38 

-44.42  (comer) 

-28.04 

•All  coefficients  in 

db . 
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It  had  been  hoped  that  the  AMP  grid  coupling  could  give 
an  indication  of  the  effects  of  a flat  plate  edge  (as  near 
the  edge  of  a wing)  on  the  coupling  coefficient.  The  greatly 
lower  values  of  coupling  coefficient  calculated  by  AMP  indi- 
cates that  the  surface  is  probably  not  being  modelled  suffi- 
ciently accurately  to  give  valid  near  field  results.  A finer 
grid  is  probably  needed  to  improve  the  quantitative  results. 

The  grid  coupling  results  show  a difference  of  28-36dB.  The 
grid  is  not  acting  as  a ground  plane. 

While  the  quantitative  results  for  this  trial  were  less 
than  encouraging,  the  qualitative  information  is  of  more 
interest.  The  results  from  equation  (10)  were  included  to 
indicate  how  it  breaks  down  calculating  near  field  information. 
The  effective  aperature  assumption  for  this  equation  makes  it 
invalid  at  these  close  distances.  The  three  positions  of 
coupling  calculation  over  the  grid  were  at  a point  equidistant 
to  an  edge  (.25m),  on  the  middle  of  an  edge  (.50m),  and  on  a 
corner  (,707m).  Using  the  "on  plate"  point  (.25m)  as  a base, 
the  edge  coupling  point  had  dropped  6.48dB  as  opposed  to  a 
change  for  that  distance  in  the  ground  plane  calculations  of 
4.72dB.  This  result,  that  coupling  was  decreased  by  moving 
the  receive  antenna  near  an  edge,  was  expected.  However,  the 
corner  data  point  shows  increased  vice  a further  decrease  in 
coupling.  Again,  using  the  ,25m  point  as  a base,  the  corner 
coupling  coefficient  drops  only  .43dB  whereas  the  same  distance 
change  over  a ground  plane  causes  a drop  of  6.84dB. 
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Again  it  is  emphasized  that  no  firm  faith  can  be  placed 
in  the  accuracy  of  these  results  in  predicting  a real  world 
situation.  The  grid  spacing  is  obviously  not  correct  to 
properly  model  a solid  plate  in  the  near  field.  The  corner 
increase  in  coupling  could  be  due  to  a corner  diffraction 
effect  or  simply  to  a built-in  anomoly  in  the  model. 

Clearly,  more  runs  need  to  be  studied  with  varied  geom- 
etries to  test  the  above  results.  As  a first  step,  diagonal 
cross  elements  in  the  vicinity  of  the  receive  antenna  or  a 
finer  grid  should  be  tried.  However,  it  is  to  be  considered 
that  these  runs  with  2S0  elements  averaged  29  minutes  of 
computer  time  on  the  NPS  IBM  360/67.  Due  to  this  time  problem, 
furthei  studies  were  not  attempted  in  this  area. 

C.  AMP  LIMITATIONS 

The  greatest  limitation  on  AMP  comes  from  the  restriction 
to  model  wire  segments  in  lengths  of  one-tenth  wavelength  or 
less  and  to  model  surfaces  with  a grid  of  the  same  order  to 
obtain  reasonably  accurate  results.  Thus,  surfaces  of  one 
square  wavelength  require  at  least  200  segments,  a one  wave- 
length cube  at  least  1200  segments,  etc.  As  can  be  seen,  the 
problem  gets  out  of  hand  quickly.  The  AMP  program  has  been 
issued  in  three  versions  allowing  either  250,  500,  or  1000 
segments  depending  on  the  accuracy  of  the  machine  on  which 
it  is  to  be  used.  There  is  a fundamental  limitation  on  the 
size  of  matrices  which  modern  computers  can  invert  without 
excessive  loss  of  accuracy.  So,  as  a rule  of  thumb,  AMP 
should  not  be  applied  to  complex  structures  of  greater  than 
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one  or  two  wavelengths  in  size.  Obviously,  applying  this 
restriction  to  the  300  megahertz  frequency  of  this  study 
limits  application  to  shapes  on  the  order  of  one  to  two 
meters  --  far  from  the  desire  of  modelling  an  aircraft. 

Figure  1 shows  some  typical  run  times  of  AMP  programs 
for  various  numbers  of  segments.  The  NPGS  computer  terminal 
time  accounting  routine  is  only  accurate  to  plus  or  minus 
about  ten  percent  depending  on  other  jobs  in  the  system. 

Several  schemes  were  investigated  to  try  to  increase 
the  size  that  AMP  could  handle.  Using  coordinate  reflection 
techniques  to  maintain  matrix  symmetry  only  yielded  an  in- 
crease of  a factor  of  four  for  a full  reflection  in  two  axes. 

An  interesting  technique  was  found  to  eliminate  unwanted 
structure  caused  by  the  reflection.  The  unwanted  portions 
can  be  connected  together  two  at  a time  with  net  work  (NT) 
cards  that  have  zero  values  entered  and  the  segments  disappear 
as  far  as  affecting  the  solution. 

Tests  were  run  with  larger  spacings  and  segment  size  to 
pet  a handle  on  the  onc-tenth  wavelength  restriction.  A two 

- v 

wavelength  square  was  modelled  with  wires  around  the  outer 
edge,  across  the  diagonals  and  connecting  the  midpoints  of 
opposite  sides.  The  wires  were  split  into  segments  of  one- 
tenth  wavelength.  The  feed  was  placed  one-quarter  wavelength 
above  the  center  of  the  "plate."  The  resulting  radiation 
pattern  indicated  that  the  plate  was  almost  totally  transparent! 
Various  other  tests  were  run  and  all  confirmed  the  one-tenth 
wavelength  limitation  for  this  plate  problem. 
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One  notable  exception  to  the  transparency  problem  of 
coarse  grid  models  was  a wing  model  that  shows  a fair  degree 
of  correlation  to  results  of  a similar  wing  from  the  BENT 
program  explained  in  the  next  section.  The  results  surprised 
this  investigator  and  are  presented  to  show  one  of  the  pitfalls 
lurking  in  this  approach.  Some  structures  can  be  modelled  in 
a very  crude  manner  and  ‘'reasonable"  results  achieved,  whereas 
other  structures  even  with  virtually  exact  modelling  will  pro- 
vide erroneous  data  as  with  the  Bell  47G-4A  helicopter  in 
[Ref.  4].  The  segmented  wing  is  shown  in  Figure  2.  The  results 
for  AMP  and  BENT  are  shown  in  Figure  3.  The  field  cut  is  in 
the  Y-Z  plane  with  zero  degrees  along  the  Z axis. 

Attempts  were  also  made  to  model  the  edge  of  a wing  using 
the  cliff  or  wire  screen  options  in  AMP.  None  of  these  ap- 
proaches proved  viable.  The  cliff  problem  only  affected  far 
field  radiation  pattern  results  and  the  wire  screen  had  to  be 
placed  on  a ground  plane. 

D.  GEOMETRICAL  THEORY  OF  DIFFRACTION  (GTD) 

W.  Dennis  Burnside  of  Ohio  State  University  has  done 
extensive  work  using  the  GTD  technique  to  study  radiation 
patterns  on  aircraft.  This  work  has  included  a series  of  Navy 
contracts  which  have  resulted  in  two  computer  programs. 

The  original  effort  was  to  produce  a computer  program  that 
would  nearly  exactly  model  an  aircraft  shape  with  cylinders, 
cones,  portions  of  spheres  and  plates.  Also,  attempts  have 
been  made  to  model  the  fuselage  with  a convex  surface  of 
revolution  [Refs.  7 and  8]. 
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BENT  IS  SOLID  LINE,  RMP  BROKEN 

fig. 3,  BENT/AMP  wing  comparison 


Results  with  shapes  other  than  plates  and  elliptic  cross 
section  infinite  length  cylinders  have  been  less  than  encour- 
aging thus  far.  Two  computer  programs  are  now  held  by  the 
Naval  Air  Development  Center,  Warminster,  Pa.  These  are  named 
ROLL  and  BENT.  ROLL  models  an  aircraft  with  an  elliptic 
cylinder  for  the  fuselage  and  plates  for  wings.  BENT  is  a 
simplification  of  ROLL,  examining  only  flat  or  bent  plate 
situations.  Each  program  had  a potential  application  in  the 
coupling  problem  and  therefore  copies  were  obtained  from  NADC. 

It  is  appropriate  to  note  here  that  there  was  no  documen- 
tation done  on  these  programs  save  some  cc  .ment  cards  in  the 
program.  Also,  the  programs  had  been  written  at  Ohio  State 
for  an  IBM  360  and  then  modified  for  use  on  the  NADC  CDC6600. 
Also,  modifications  had  been  made  at  NADC  for  different  output. 

i 
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| The  NPS  computer  is  an  IBM  360/67  but  the  tape  from  NADC  was 
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in  a different  character  set,  had  incompatible  typed  output 
and  pen  plotter  output.  Without  documentation,  the  cleaning 
up  of  these  programs  for  use  on  the  NPS  computer  proved  tedious 
and  time  consuming.  Appendix  A is  included  to  shed  some  light 
on  use  of  these  programs  for  future  users. 

As  indicated  in  the  introduction,  these  programs  use  a ray 
optics  technique  to  solve  for  a field.  The  primary  components 
of  the  energy  at  a point  in  a field  are  direct,  reflected  and 
diffracted  energy.  The  approach  in  obtaining  the  field  for  a 
• structure  is  to  calculate  all  the  components  (direct,  reflected 

f and  diffracted)  from  each  surface  at  one  point  in  space  and  sum 
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The  frequency  range  of  application  of  this  method  is  not 
clear.  Input  is  given  in  gigahertz.  Reference  [7]  calls  it 
a high  frequency  technique  and  [Ref.  8]  includes  information 
on  a light  pen  adaptation  of  the  program  with  selectable 
frequencies  from  one  to  ten  gigahertz.  Discussion  with  various 
individuals  indicated  that  the  GTD  approach  was  limited  to 
microwave  frequencies. 

It  was  originally  hoped  that  the  programs  could  be  modified 
to  give  near  field  and  near  structure  field  strengths.  This 
has  not  been  accomplished  to  date.  One  initial  problem  in  this 
area  is  that  the  source  and  various  scattering  centers  be  sep- 
arated by  at  least  a wavelength.  This  restriction  ties  in  with 
the  lowest  usable  frequency  discussion  in  the  previous  para- 
graph. 

The  sample  programs  in  Appendix  A give  some  indication  of 
the  input/output  combinations  available. 

It  is  interesting  to  compare  the  plot  outputs  of  Appendix 
A.  The  ROLL  solution  includes  the  effects  of  the  infinite 
cylinder  whereas  BENT  does  not.  Past  forty-five  degrees  from 
vertical,  the  patterns  are  virtually  identical.  The  sharp 
step  at  about  100  and  260  degrees  is  probably  caused  by  the 
program  switching  calculating  schemes  from  direct  and  reflected 
to  diffracted.  These  steps  seem  to  occur  more  often  when  there 
are  shallow  grazing  angles  from  source  to  edge.  This  indicates 
that  the  diffraction  coefficients  used  are  not  precise  at 
shallow  angles  because  they  include  only  first  order  effects. 
Generally,  the  true  field  pattern  appears  to  be  a faired  curve 
through  the  jump. 
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E.  BENT/ AMP  COMPARISONS 

As  discussed  earlier,  BENT  is  supposedly  restricted  from 
use  at  other  than  microwave  frequencies  and  with  a feed  closer 
than  one  wavelength.  This  limitation  is  due  to  the  inaccuracy 
of  the  diffraction  coefficients  at  low  angles.  However,  these  1 

I 

restrictions  were  discovered  after  successfully  running  a two  ! 

meter  flat  plate  with  a feed  at  .1  meter  from  the  center  at 
300  megahertz.  It  seemed  to  give  results  that  were  reasonable. 

As  discussed  in  section  II.  C.  above,  it  was  shown  that 
for  reasonable  results  with  AMP,  the  surface  must  be  modelled 
with  a grid  that  has  segments  and  holes  no  larger  than  .1  wave- 
length. This  means  that  there  must  be  about  200  segments  plus  j 

i 

final  edge  segments  per  square  wavelength  of  surface.  Also,  I 

as  the  segments  in  AMP  increased  over  200  the  run  time  became  j 

5 

20  to  30  minutes  and  two -day  turnarounds  became  normal  at  the  ; 

i 

NPS  computer  facility.  Thus  at  300  MHZ  the  maximum  realistic  j 

surface  was  one  square  meter. 

It  was  decided  to  run  a comparative  study  of  the  BENT  and 

AMP  programs  to  get  a better  grasp  on  the  above  limitations. 

A one-square  meter  plate  was  modelled  with  both  BENT  and  AMP. 

The  feed  was  placed  at  .08,  .5,  and  1.0  meter  positions  above 

the  center  of  the  plate  as  in  Figure  4.  In  Figures  5 through  j 

8 the  results  of  this  comparison  are  shown. 

Reasonably  high  confidence  is  placed  in  the  AMP  results 

due  to  the  fineness  of  the  grid  and  the  past  history  of  AMP  j 

data.  The  feed  in  AMP  was  a .1  meter  whip  divided  into  ten  j 

segments  and  fed  from  the  lower  segment.  It  was  positioned  \ 
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Geometry  of  Flat  Plate 
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E-PHI  DB  PLOT  40  DB  SC, RLE 

NO.  OF  CORNERS  OF  PL  RTF.--  4 

ROTRTION  RNGLF  GF  PLRTE---  180.000 

FREQUENCE  0.300 

RNTENNR  TVDE---  0.08N  FEED 

THETR  ROTRTION  RNGLE  OF  DIPOLE---  90-000 

PHI  ROTRTION  RNGLE  OF  DIPOLE-  0.000 

PATTERN  NO.  PLRTE 

BENT  IS  SOLID  LINE,  RMD  BROKEN 

fig. 5,  Flat  plate  BENT/AMP  comparison,  ,08m  feed 


E-PH]  DB  PLOT  40  DB  SCALE 
NO.  OF  CORNERS  OF  PlRTf>  4 
ROTRTION  RNGLE  OF  PLRTE--  150-000 
FREQUENCE-  0-300 
RNTENNR  T^PE=  . 5M  FEED 

THETR  ROTRTION  RNGLE  OF  DIPOLE--  90-000 
PM]  ROTRTION  RNGLE  OF  DIPOLE----  0-000 
PRTTERN  NO.  PLATE 
BENT  IS  SOL  I D LINE,  RMP  BROKEN 

fig. 6.  Flat  plate  BENT /AMP  comparison,  .5m  feed 
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so  that  the  point  .03  meter  from  the  bottom  was  at  .08,  . 5 or 
1.0  meter.  The  BENT  source  was  an  infinitesimal  dipole. 

The  jumps  in  the  BENT  solutions  usually  occur  when  the 
program  shifts  methods  of  calculating  the  field.  A good  example 
of  this  is  the  step  at  1S3.S  and  206.5  degrees  on  the  .S  meter 
feed  plot.  Geometry  of  the  problem  shows  that  this  is  exactly 
where  the  region  shadowed  from  the  feed  begins.  Usually  these 
bumps  in  the  BENT  solution  can  be  just  faired  through.  It  is 
important,  however,  to  be  sure  that  the  bump  is  from  extraneous 
causes  and  not  real. 

The  comparative  results  show  that  in  this  case  the  two 
methods  are  in  genial  agreement.  With  the  exception  of  the 
lower  lobe  nuances  at  .08  m feed  and  the  magnitude  of  the  upper 
lobes  at  . 5 m feed,  the  two  methods  show  surprising  agreement. 

It  is  interesting  to  note  also  that  the  average  computer 
run  time  for  BENT  including  compiling  was  less  than  two  min- 
utes (about  50  seconds  compile  time)  and  the  average  run  time 
of  AMP  which  is  precompiled  was  twenty-seven  minutes. 

The  results  indicate  that  BENT  will  give  usable  data  at 
UHF  frequencies  and  at  locations  within  a half  wavelength  of 
the  surface.  For  small,  square  plates,  apparently  BENT  may 
be  applied  with  feeds  essentially  on  the  surface  (.1  wave- 
length) . 

F.  TCI  RESULTS 

Robert  Tanner  of  Technology  for  Communications  International 
was  kind  enough  to  calculate  the  one-meter  plate  problems  of 
section  II.  E.,  using  the  TCI  wire  grid  solving  computer  program. 
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His  solutions  are  presented  over  the  BENT/AMP  solutions  in 
Figures  8 through  10. 

Tanner  expressed  great  confidence  in  these  results  due  to 
past  experience  with  the  TCI  program. 

The  exact  formulation  of  the  TCI  program  is  confidential 
within  TCI  but  it  is  known  that  a wire  grid  technique  is  used. 
However,  the  grid  spacing  can  be  much  larger  than  .1  wavelength. 
This  is  accomplished  by  recognizing  that  if  a wire  square  is 
used  to  model  a square  patch  of  surface  that  the  wire  'loop' 
will  have  more  inductance  and  less  capacitance  than  the  plate. 
The  TCI  program  takes  this  effect  into  account  whereas  AMP 
does  not.  Run  times  on  the  TCI  CDC6600  computer  were  about 
5 seconds. 

G.  HYBRID  TECHNIQUES 

Reference  9 shows  one  method  of  combining  MM  and  GTD  to 
expand  these  methods  to  problems,  including  a monopole  next  to 
a conducting  edge  and  a monopole  on  a plate.  The  approach  was 
to  start  with  a MM  formulation  and  modify  the  impedance  matrix 
based  on  GTD  considerations.  The  results  show  that  the  input 
impedance  and  reactance  of  a monopole  near  a conducting  edge 
of  90  degrees  (cliff)  begin  to  oscillate  above  and  below  the 
input  values  on  an  infinite  ground  plane  when  the  monopole  is 
within  about  three  wavelengths  of  the  edge.  The  oscillations 
become  most  pronounced  inside  one  wavelength  and  the  approach 
generally  broke  down  inside  .2  wavelengths.  Similar  results 
were  chieved  for  a monopolc  over  an  octagon,  square  and 
circular  plate.  The  octagon  results  arc  better  than  the 
square  since  vertex  diffraction  is  neglected. 
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FIG. 8.  BENT/AflP  PLATE  WITH  TCI  RESULTS.  .08m  FEED 
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fig. 9.  BENT/AMP  plate  with  TCI  results,,  .5m  feed 
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E-PHI  OB  PLOT  40  DB  SCALE 
HO.  OF  CORNERS  OF  PLATE--  4 
ROTATION  ANGLE  OF  PLATE-  180.000 
FREQUENCE  0-300 
ANTENNA  TVPE--  1H  FEED 

THETA  ROTATION  ANGLE  OF  DJP0LE-  90-000 
' PHI  ROTATION  ANGLE  OF  DIPOLE-  0-000 

PATTERN  NO.  PLATE 

BENT  IS  SOLID  LINE,  AMP  BROKEN 7TJ  ASTER JCK 

FIG. 10.  BENT/AMP  PLATE  WITH  TCI  RESULTS/  I.Om  FEED 
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it  was  hoped  to  obtain  far  field  data  or*  this  work  to 
compare  with  tne  plate  calculations  in  section  II.  E.  and  to 
attempt  to  calculate  coupling  coefficients  with  this  technique. 
The  conducting  wedge  and  plate  approaches  seem  ideal  to  model 
wing  surfaces.  However,  due  to  time  limitations  this  work  has 
not  been  accomplished. 
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III.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  solution  of  antenna  performance  near  complex  conducting 
structures  is  still  a problem.  Computer  techniques  using  moment 
methods,  geometrical  theory  of  diffraction  and  modifications 
and  combinations  thereof  are  useful  for  certain  specific  and 
generally  simple  problems,  but  no  universal  solution  method 
is  now  available. 

The  MM  technique  is  too  cumbersome  for  application  to 
structures  much  larger  than  one  wavelength  and  the  OSU  GTD 
programs  cannot  be  applied  to  the  coupling  problem  since  there 
are  no  near-field  data  (field  within  the  size  of  the  structure). 

The  data  obtained  from  ground  plane  analysis  using  AMP 
indicates  1.24  db  less  coupling  for  real  world  UHF  antennas 
than  theory  predicts.  These  results  may  be  helpful  in  design- 
ing antenna  arrays  on  large  flat  surfaces.  However,  these 
results  have  not  been  substantiated  with  actual  measurements 
so  they  may  still  be  optimistic  or  pessimistic. 

The  most  promising  method  for  solution  of  coupling  problems 
appears  to  be  a hybrid  approach  using  GTD  to  modify  the  im- 
pedance matrix  in  an  MM  formulation.  TCI’s  "modified”  surface 
approach  modifies  the  matrix  in  a similar  manner.  This  method 
readily  gives  the  current  in  the  various  segments  and  allows 
easy  calculation  of  the  coupling  coefficients. 

However,  even  with  the  most  recently  published  hybrid 
methods,  the  solution  of  edge  fields  has  yet  to  be  addressed. 
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Also,  more  exact  diffraction  coefficients  need  to  be  developed 
to  allow  closer  placement  of  sources  and  diffraction  surfaces. 
The  included  results  for  a square  plate  show  that  in  some 
cases,  however,  satisfactory  results  may  be  achieved  with 
present  coefficients.  As  the  various  methods  become  more 
refined,  higher  order  terms  that  have  so  far  been  neglected 
will  have  to  be  considered,  such  as  double  diffraction  and 
sharp  angle  corner  diffraction. 
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APPENDIX  - OSU  GTD  PROGRAMS  "ROLL"  AND  "BENT" 


I.  ROLL  PROGRAM 

A.  GENERAL 

The  OSU  ROLL  program  combines  an  infinite  cylinder  of 
elliptic  cross  section  with  "wings"  which  may  be  bent  once. 

The  fields  are  calculated  from  contributions  due  to  direct 
source  energy,  reflected  energy  and  diffracted  energy.  The 
program  calculates  the  fields  from  one  wing,  then  adds  on  the 
fields  from  succeeding  wings. 

B.  COMMENTS  ON  INPUT  DATA 

See  Table  A-I  for  a listing  of  input  data  and  Figure  A-l 
for  the  coordinate  system  used. 

It  was  found  that  THC  and  PHC  cards  had  to  be  entered  at 
some  value  other  than  zero,  at  least  with  the  geometry  of  the 
sample  problem.  If  entered  as  zero,  the  program  would  fail  in 
subroutine  SGTD  on  card  011450  in  calculating  TAN  (AS).  AS 
was  90  degrees.  At  this  point  SGTD  was  branched  from  FUSLAG 
which  was  branched  from  REFFD  which  was  called  from  MAIN. 

Note  that  wing  geometry  must  not  penetrate  fuselage.  It  is 
important  that  geometries  closely  coincide.  They  should  be 
checked  on  a calculator.  The  NPGS  computer  is  an  IBM  360. 

C.  OUTPUT 

The  version  of  the  program  received  at  NPS  had  various 
minor  output  incompatabilities.  Notably,  all  quotations 
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were  double  and  had  to  be  changed  to  single.  Also,  some 
foresting  of  data  was  awkward,  such  as  titles  following  im- 
mediately below  data  and  then  a page  skipped  prior  to  actual 
data  listing.  The  POLPLT  subroutine  had  to  be  rewritten  to 
make  it  compatable  with  the  NPS  plot  subroutine  package. 
Fortunately,  the  ROLL  authors  and  NPS  both  generally  followed 
Caicomp's  recommendations  in  plot  subroutines. 

The  POLPLT  subroutine  plots  either  magnitude  power  or  db 
polar  plots.  The  plots  are  scaled  to  fit  all  data  into  the 
plot  size  indicated  by  RADC.  In  the  main  part  of  the  program, 
E-theta  and  E-phi  values  are  searched  for  the  maximum.  This 
is  then  called  EMX.  Within  the  POLPLT  subroutine,  RADC  is 
called  RP  and  E-theta  or  E-phi,  as  appropriate,  is  called  ET 
(which  is  complex).  An  intermediate  terra  is  formed,  ETM, 
which  is  the  complex  absolute  value  of  ET  divided  by  EMX.  The 
radius  of  the  point  on  the  polar  plot  RD  is  found  as  follows: 
For  Field  Plot: 

RD  “ RP  x ETM 
For  Power  Plot: 

RD  * RP  x ETM  x ETM 
For  db  Plot: 

RD  * 2D  x A LOG 10  (ETM) 

IF  (P.D.LT.  - 40)  RD  * -40 
RD  - RP  x (RD  ♦ <0)/40 

The  logical  LWRT  (card  3)  controls  printing  of  intermediate 
values  of  the  E-theta  a ltd  E-pL\  field.  Thus,  for  one  wing, 
this  is  722  lines  of  output  --  for  two  wings,  1444,  etc.  The 

final  listings  of  the  field  are  not  affected  by  LWRT. 
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Card  5 allows  skewing  of  the  field  pattern  axis.  Note 
that  ITHI,  ITHF,  etc.,  apply  in  this  "new”  coordinate  rotation. 

D.  GEOMETRY 

If  the  wings  are  located  below  the  centerline  of  the 
fuselage,  a tangent  fillet  from  the  widest  point  in  the  fu- 
selage is  dropped  down  to  the  wing. 

The  program  "blows  up"  if  the  wings  are  placed  on  top  of 
the  fuselage  with  the  source  also  on  top,  according  to  NADC. 

E.  SAMPLE  DATA  DECK 

The  data  deck  in  Figure  A-2  is  fcr  a mid-wing  model  with 
the  fuselage  72  inches  in  diameter.  The  overall  wing  span  is 
twenty  meters  with  a three-meter  root  and  even  tapering  to  a 
one-meter  tip.  The  wing  is  flat. 

The  source  is  the  monopole  and  is  located  exactly  over  the 
center  of  the  wing  on  the  top  of  the  fuselage.  The  monopole 
is  oriented  perpendicular  to  the  plane  of  the  wings.  See 
Figure  A-3  for  a sketch  of  the  sample  problem.  The  field 
plot  data  is  taken  in  the  X-Y  plane  with  the  zero  degree 
position  corresponding  to  the  positive  X axis. 

The  polar  plot  output  is  shown  in  Figure  A-4. 


TABLE  A- I 

Input  Format  for  ROLL 


Card  Column 
Number  Number 


6-10 

11-15 


11-15 

16-20 


Description 


MP  ■ Number  of  wings 

MX  • Number  of  edges  per  wing 

FRQG  ■ Frequency  in  gigahertz 


AF  - One-half  overall  X-direction 
of  fuselage 

BF  * One-half  overall  Y-direction 
of  fuselage 

PHSO  ■ Angular  position  of  source, 
in  degrees,  zero  equals  top 
of  fuselage 

ZS  = Z-coordinate  of  source 
JANT  * Type  of  antenna 

(infinitesimal  element) 
slot  ■ 1 

radial  monopole  2 
BETA  = Angle  slot  makes  with  Z-ax:s 

Note:  Distances  entered  in  inches 


LWRT  ■ Additional  write-out  desired 
(true  or  false) 

LPLT  ■ Pattern  plotted  by  pen  plotter 
(true  or  false) 

LFLAT  « Is  flat  plate  considered 
(true  or  false) 


THETA  and  PHI  are  radiation  angles 
ITHI  ■ Initial  value  of  THETA  in 
degrees 

5-10  I ITHF  * Final  value  of  THETA  in 
degrees 

ITHS  * Incremental  step  in  THETA  in 
degrees 

IPHS  c Incremental  step  in  PHI  in 
degrees 

Note:  If  ITHI  does  not  equal  ITHF 

then  cards  after  card  5 must 
be  reentered  for  eac''!  step  in 
THETA.  Also,  as  the  program 
is  set  up,  IPHS  * 1 


Format 


2110, 

F10. 5 


4F10.5 

110, 

F10.5 


ITHS 


Note: 


TABLE  A-l  (cont'd) 


Card 

Number 


Column 

Number 


1-10 

li-20 


1-10 

11-20 

21-30 


7 

to 

6+MX 


1-10 

11-20 

21-30 


7+MX 


1-10 

11-12 


Description 

Format 

THC  = THETA  rotation  angle  of  the 
axis  of  pattern  rotation 
PHC  = PHI  rotation  angle  of  the 
pattern  rotation 

2F10.S, 

Note:  Angles  in  degrees.  ITHI 

and  ITHF  are  values  in  this 
rotated  coordinate  system 

MCI  = Number  of  first  corner  at  bend 
MC2  = Number  of  second  corner  at  bend 
RTH  = Angle  of  bent  between  surfaces 

2110 

F10.5 

Note:  If  wing  is  defined  flat 

(LFLAT=T) , then  MC1=MC2=MX 
. and  RTH=180 

X coordinate  of  a corner  of  the  wing 
Y coordinate 
Z coordinate 

3F10.S, 

Note:  Positions  entered  in  inches. 

Points  must  be  ordered  counter 
clockwise,  first  point  on 
fuselage.  At  no  point  may  the 
wing  plane  enter  the  fuselage 

Only  used  if  LPLT  = T 
RADC  = Radius  of  polar  plot  on  pen 
plotter 

I PLT  = 

Field  Plot  = 1 
Power  Plot  = 2 
DB  Plot  = 3 

rlO.5, 

12 

If  more  than  one  wing  is  called  for 
(MP  greater  than  1),  then  repeat 
cards  6 and  7 to  MX  with  the  appro- 
priate coordinates  for  new  wing. 

If  ITHI  does  not  equal  ITHF,  then 
all  cards  from  first  card  6 on 
| must  be  repeated  for  each  step  in 
j THETA. 

123456789  FOLLOWING  IS  SAMPLE  ROLL  INPUT  DECK 


36.  36. 

TTTTTTTTTTTTTTT 


.01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

3. 


90 


90 


.01 

♦ 

-36. 
-393.7 
-393.7 
-36  • 

* 

36. 

393.7 

393.7 

36. 

3 


4.3 

0. 

1 

4180. 

59. 

19.7 

-19.7 

-59. 

4180. 

-59. 

-19.7 

19.7 

59. 


0. 


123456789  FOLLOWING  IS  SAMPLE  BENT  INPUT  DEC< 

TTTTTTTTTTTTTTTWING  FLAT  36INCH  FEED 
0 2 70  90 


.01 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

1. 

36. 

90. 

3. 


90 


90  1 

1.0 
6.3 
0. 

-393.7 

-393.7 

0. 

393.7 

393.7 

0. 

0. 

0. 

0. 

3 


1 


180. 

59. 

19.7 

-19.7 

-59. 

-19.7 

19.7 


0. 


fig.A-2.  Sample  data  decks  for  ROLL  and  BENT 
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II.  BENT  PROGRAM 


A.  GENERAL 

The  OSU  BENT  program  calculates  the  far  field  from  z flat 
plate  with  one  bend  in  it.  The  fields  are  calculated  from 
contribution  due  to  direct  source  energy,  reflected  energy 
and  diffracted  energy.  The  program  is  a simplification  of 
the  ROLL  program  only  considering  one  "wing.”  However, 
multiple  sources  (NADC  said  up  to  three  but  the  program  seems 
to  allow  any  number --not  checked)  may  be  handled.  The  pro- 
gram has  application  when  modelling  flat  fuselage  bottoms, 
wings  or  tail  areas. 

The  READ  (5,9966)  RADC , IPLT  card  (#005120)  had  to  be 
separated  from  the  call  plots  card  to  run  on  the  NPS  IBM  360. 

B.  COMMENTS  ON  DATA  INPUT 

Table  A-II  lists  detailed  card  entries  for  BENT. 

For  the  orientation  in  the  sample  problem,  ANGS  had  to  be 
entered  at  other  than  zero  so  that  the  final  value  of  E-phi 
would  not  be  excessive  and  mis-scale  the  pen  plots.  NADC 
says  that  RTH  equal  to  90  will  not  work.  This  has  not  been 
checked.  Small  (one  lambda  or  less)  plates  seem  to  take 
excessive  run  time  and  produce  a size  warning  in  the  program. 

C.  OUTPUT 

The  POLPLT  subroutine  required  extensive  rewriting  for 
comparability  with  the  NPS  plot  routines.  Other  output  and 
control  are  similar  to  ROLL. 
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D.  GEOMETRY 


The  ci  rdinate  system  is  as  used  in  ROLL,  as  shown  in 
Figure  A-l. 

E.  SAMPLE  DATA  DECK 

The  sample  problem  for  BENT  is  just  the  wing  of  the  ROLL 
sample  problem  without  the  fuselage.  The  wing  is  twenty  meters 
in  span,  three  meteis  at  the  root,  and  one  meter  at  the  tips. 
The  feed  is  a dipole  directly  over  the  center  cf  the  wing 
thirty-six  inches  (.914m)  from  the  wing  and  oriented  perpen- 
dicular to  the  plane  of  the  wing.  In  other  words,  the  feed 
is  in  exactly  the  same  location  as  on  ROLL  except  the  fuselage 
is  missing.  The  sample  data  deck  is  listed  in  Figure  A-3  and 
polar  plot  data  similar  to  that  for  ROLL  is  in  Figure  A-5. 
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TABLE  A- I I 

Input  Format  for  BENT 

Card 

Column 

Number 

Number 

Description 

Format 

1 

1-5 

LWRT 

= Additional  write-out  desired 

3L5 , 

(true  or  false) 

6-10 

LPLT 

= Plot  pattern  with  pen  plotter 

2A5 , 

(true  or  false) 

11-15 

Ll'LAT 

= Is  flat  plate  considered 

A10 

(true  or  false) 

16-20 

CRAFT 

= Label  for  plot,  aircraft  type 

21-25 

PATT 

= Label  for  plot,  pattern  number 

26-35 

TYPE 

= Label  for  plot,  antenna  type 

2 

1-5 

PT  = 

Label  top  of  plot 

3A5 

6-10 

PL  = 

Label  left  of  plot 

11-16 

PR  = 

Label  right  of  plot 

THETA 

and  PHI  are  radiation  angles 

i 

5 

1-5 

ITHI 

= Initial  value  of  THETA  in 

515 

degrees 

6-10 

I THE 

= Final  value  of  THETA  in 

degrees 

11-15 

ITHS 

= Incremental  step  in  THETA 

I PHS 

= Incremental  step  in  PHI 

NUM 

= Number  of  sources 

Note : 

If  ITHI  does  not  equal  ITHF, 
then  the  three  source  cards 

(per  source)  and  the  plot  size/ 
type  card  must  be  reentered 
for  each  step. 

As  written,  IPHS  must  equal  1 

4 

1-10 

ANGS 

= Initial  value  of  PHI 

2F1 0 . 5 

11-20 

DELPHI  = Incremental  value  of  PHI 

Note : 

As  written,  DELPHI  must  equal 
1.0 

5 

1-5 

MX  = 

1 

Number  of  corners 

215, 

G-10 

MC  = 

Number  of  corner  junctions, 
with  bend  between  corner  H 1 

2F1 0 . 5 

and  MC 

13-20 

FRQG 

= Frequency  in  gigahertz 

21-30 

RTH 

= Angle  of  bend  in  plate,  only 

defined  between  90  and  270  degrees 

Note  : 

If  LFLAT  is  true,  then  MC=MX 
and  RTH  = 180.  See  input  discussion 
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TABLE  A- I I (cont’d) 


Card 

Number 

Column 

Number 

Description 

Format 

6 

to 

5+MX 

HI 

Corner  locations  in  inches 
X coordinate 
Y coordinate 
Z coordinate 

3F10.S 

i 

6+MX 

1-10 

THC  ■ THETA  rotation  angle  of  the 
axis  of  pattern  rotation 
PHC  “ PHI  rotation  angle  of  the 
axis  of  pattern  rotation 

2F10.5 

Note:  Angles  in  degrees.  ITH1 

and  ITHF  are  values  in  this 
rotated  coordinate  system 

The  next  three  cards  control  the 
location  and  orientation  of  sources. 
Repeat  these  cards  for  each  source. 

'/♦MX 

1-10 

11-20 

ACUR  K Amplitude  of  source 
PHCUR  = Phase  of  source  in  degrees 

2F1U.5 

8+MX 

rai 

II  Erl™ 

HWfi» 

Location  of  source  in  inches 
X coordinate 
Y coordinate 
Z coordinate 

3F10.5 

9+MX 

1-10 

11-20 

THOO  = THETA  rotation  angle  of 

source  with  respect  to  ref- 
erence coordinate  system 
PHOO  = PHI  relation  angle  of  source 

2F10.5 

i n j iiv 
A l/T  J'lA 

1-10 

11-12 

Use  only  if  LPLT  is  true 
RAI3C  - Radius  of  polar  plot  on  pen 
plotter  in  inches. 

Recommend  3 

IFLT  ■= 

1 ■*  Field  Plot 

2 = Power  Plot 

3 - DB  Plot 

F10. 5 , 
12 

: 

! 

Note:  Again,  if  there  are  steps  in 

THETA,  then  the  source  and 
plot  cards  must  be  repeated 
for  each  step. 

WING 


0 


E-PHJ  DB  PLOT  40  DB  SCALE 

HO.  OF  CORNERS  OF  PLATE---  6 

ROTATION  ANGLE  OF  PLATE---  180.000 

FREQUENCE-  0*300 

RNTENNR  TVPE~  38INCH  FEE 

THETR  ROTATION  ANGLE  OF  DIPOLE--  90-000 

PHI  ROTATION  ANGLE  OF  DIPOLE-  0-000 

PATTERN  NO-  FLAT 

fig.A-5,  BENT  output,  sample  problem 
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